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ABSTRACT 

The radially averaged metallicity distribution of the ISM and the young stellar population of 
a sample of 20 disk galaxies is investigated by means of an analytical chemical evolution 
model which assumes constant ratios of galactic wind mass loss and accretion mass gain to 
star formation rate. Based on this model the observed metallicities and their gradients can be 
described surprisingly well by the radially averaged distribution of the ratio of stellar mass to 
ISM gas mass. The comparison between observed and model predicted metallicity is used to 
constrain the rate of mass loss through galactic wind and accretion gain in units of the star 
formation rate. Three groups of galaxies are found: galaxies with either mostly winds and 
only weak accretion, or mostly accretion and only weak winds, and galaxies where winds are 
roughly balanced by accretion. The three groups are distinct in the properties of their gas disks. 
Galaxies with approximately equal rates of mass-loss and accretion gain have low metallicity, 
atomic hydrogen dominated gas disks with a flat spatial profile. The other two groups have 
gas disks dominated by molecular hydrogen out to 0.5 to 0.7 isophotal radii and show a radial 
exponential decline, which is on average steeper for the galaxies with small accretion rates. 
The rates of accretion (< 1.0 x SFR) and outflow (< 2.4 x SFR) are relatively low. The latter 
depend on the calibration of the zero point of the metallicity determination from the use of 
Hu region strong emission lines. 
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1 INTRODUCTION 

The chemical composition of stars and the interstellar medium 
(ISM) reflect the evolutionary history of a galaxy from its for¬ 
mation to the present stage. While the physical processes leading 
to chemical evolution and enrichment with heavy elements such 
as star formation, nucleosynthesis, stellar winds, supernova explo¬ 
sions, galactic winds and accretion are complicated and coupled in 
a complex way, an intriguingly simple relationship between galac¬ 
tic stellar mass and the gas-phase oxygen abundance of the ISM 
or the metallicity of the young stellar population of star form¬ 
ing galaxies has been observed, the mass-metallicity relationship 
(MZR; Lequeux et al. 1979; Tremonti et al. 2004; Kudritzki et al. 
2012). This indicates that there might be a straightforward sim¬ 
plified way to describe the key aspects of chemical evolution. For 
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instance, Lequeux et al. (1979) already argued that a simple closed- 
box chemical evolution model, where galaxies are not affected by 
mass accretion or mass loss and which relates metallicity to the 
ratio of stellar to ISM gas mass, could explain the MZR of their 
sample of irregular galaxies well. Since then it has become clear 
that in a ACDM dominated universe processes of galaxy merg¬ 
ing and accretion from the intergalactic medium of the cosmic 
web and galactic winds may play a crucial role in galaxy forma¬ 
tion and evolution. Consequently, the situation is far more complex 
than in a simple closed-box model (see, for instance, Dave et al. 
201 la,b, 2012; Yates et al. 2012; Dayal et al. 2013). Moreover, 
comprehensive spectroscopic surveys from low to high redshifts 
(see Zahid et al. 2014 and references therein) clearly indicate that 
there is an empirical upper limit in galaxy metallicity for the most 
massive galaxies (Zahid et al. 2013). This is in contradiction to the 
classical closed-box model, for which the metallicity increases log¬ 
arithmically with increasing ratio of stellar mass to gas mass. On 
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the other hand, as has been shown recently by Zahid et al. (2014), 
when allowing for galactic winds and accretion the observed MZR 
can still be explained by a model where the observed metallic- 
ity is a simple analytical function of the stellar to gas mass ratio. 
Ascasibar et al. (2014) find a very similar result. They develop a 
new approach, in which even a closed box model would eventually 
reach a maximum metallicity, equal to the oxygen abundance of the 
supernova ejecta. They also point out that their approach should ap¬ 
ply to spatially resolved observations. 

The study of global galaxy metallicity as a function of global 
stellar and gas mass is only a first step to constrain galaxy evo¬ 
lution. An important extension of these investigations, as already 
noted by Ascasibar et al. (2014), is to use the information of the 
observed resolved spatial distribution of metallicity in the disks of 
star forming galaxies. As is known since a long time (e.g., Searle 
1971; Garnett & Shields 1987; Garnett et al. 1997; Skillman 1998; 
Vila-Costas & Edmunds 1992; Zaritsky et al. 1994), disk galaxies 
show negative abundance gradients with increased metallicity to¬ 
wards the centres and lower metallicity in the outskirts. Since then 
a considerable effort has been made to model the spatially re¬ 
solved chemical evolution of galactic disks (e.g., Molla et al. 1997; 
Prantzos & Boissier 2000; Chiappini et al. 2001; Fu et al. 2009; 
Pilkington et al. 2012; Mott et al. 2013). 

This is the starting point for the study presented here. With 
the indication that the global metallicity of galaxies can be re¬ 
lated to the global ratio of stellar to gas mass by relatively sim¬ 
ple chemical evolution models (Zahid et al. 2014; Ascasibar et al. 
2014; Yabe et al. 2015), we investigate whether the spatially re¬ 
solved radial distribution of metals in disk galaxies can be ex¬ 
plained in a similar way. We use the radial distributions of stel¬ 
lar mass and ISM Hi and H 2 gas mass column densities together 
with the observed metallicities and metallicity gradients of 20 well 
observed disk galaxies to find out whether the analytical chemical 
evolution models introduced below can explain the observed radial 
metallicity distribution, and whether the influence of galactic winds 
and accretion can be constrained in this way. 

This work has been triggered by the work by Ho et al. (2015) 
who investigated IFU observations of 49 local field star-forming 
galaxies and found a narrow Gaussian distribution of metallic¬ 
ity gradients characterised by a mean ± standard deviation of 
—0.39 ± 0.18 dex (R 25 is the B-band iso-photal radius). 
Ho et al. (2015) compared this observed distribution of metallicity 
gradients with that predicted by chemical evolution models from 
a set of stellar and gas mass profiles of 13 nearby galaxies. The 
latter were taken as benchmark galaxies representing local galax¬ 
ies. They concluded from the comparison of the two statistical dis¬ 
tributions that the effects of galactic winds and accretion must be 
very small. Ho et al. (2015) used a statistical approach to compare 
the observed distribution of metallicity gradients of their sample 
of IFU-studied galaxies with a model-generated sample based on 
the observed mass profiles of a second set of galaxies. We will ap¬ 
ply a different approach and will investigate one single sample of 
galaxies with a homogeneous set of metallicity and stellar and gas 
mass observations. In addition, we will not only use the metallicity 
gradients but also the absolute level of metallicity to constrain the 
effects of winds and infall for each galaxy in detail. 

We describe the sample of local star-forming galaxies and 
their observations in Section 2. The chemical evolution model is 
introduced in Section 3 and then applied in a first step in Section 4 
to the Milky Way disk, where we use a fit of metallicity and metal¬ 
licity gradient observed with Cepheids and B-stars to empirically 
constrain the metallicity yield. The next step is then to apply our 


models for a fit of the observed gas-phase ISM oxygen abundances 
of our star-forming galaxies. This is done in Section 5 and the re¬ 
sults are discussed and summarised in Section 6. 


2 GALAXY SAMPLE AND OBSERVATIONS 

The sample investigated in this work consists of 19 nearby spi¬ 
ral galaxies selected from the study by Schruba et al. (2011) and 
Leroy et al. (2008), who published carefully measured radial pro¬ 
files of Hi and H 2 gas mass column densities. We combine 
these observations with measurements of oxygen abundances and 
abundance gradients published in the comprehensive study by 
Pilyugin et al. (2014). The study by Schruba et al. (2011) contains 
33 spirals, however, of those only 21 overlap with Pilyugin et al. 
(2014). Of those 21, NGC2841 and NGC4725 were not included 
into our final sample, because the oxygen abundance observations 
covered only a small radial range of the galaxies and were also 
very uncertain. The basic properties of our sample are summarised 
in Table 1. We note that while the majority of our sample con¬ 
sists of isolated field galaxies, NGC5194 (M51) is a merger and 
NGC4625 is a pair with NGC4618 (the latter is not included in 
our sample). For further information see Schruba et al. (2011) and 
Leroy et al. (2008). In the following, we use the term “SL-sample” 
for this sample of 19 spiral galaxies. 

We include the Milky Way as the 20th galaxy in our investi¬ 
gation. We will use the Milky way as test galaxy for the chemical 
evolution models and for the empirical calibration of our yields. 

2.1 HI and H2 mass column density profiles 

2.1.1 The SL-sample 

The radial column densities of atomic hydrogen, Ehi, were ob¬ 
tained from VLA observations of the 21 cm hydrogen line car¬ 
ried out in the The Hi Nearby Galaxy Survey survey (THINGS; 
Walter et al. 2008). A factor 1.36 was included to account for he¬ 
lium and heavier elements (see Schruba et al. 2011). The molecu¬ 
lar gas column densities, E/r 2 , result from CO(2 —> 1) line emis¬ 
sion observations with the IRAM 30m telescope obtained as part 
of the HERA CO Line Extragalactic Survey survey (HERACLES; 
Leroy et al. 2009). For all details, such as the conversion from CO 
to H 2 , the azimuthal averaging and construction of deprojected ra¬ 
dial profiles, the conversion of emission line intensities into mass 
column densities, etc., the readers are referred to Schruba et al. 
(2011), who present the radial profiles of Hi and H 2 gas mass col¬ 
umn densities. The total ISM gas mass column density is then E 9 
= E h \ + S h 2 - 

2.1.2 The Milky Way 

For gas mass column densities of the Milky Way disk, we use the fit 
formulae provided by Wolfire et al. (2003). We restrict ourselves to 
the range of galactocentric distances between 4.5 to 9 kpc for which 
a good description of the stellar mass profile and accurate measure¬ 
ments of the metallicity and metallicity gradient are available (see 
below). In this range, Wolfire et al. (2003) suggest a constant mass 
column density of the atomic gas, E h\ = 6.8 Mg pc~ 2 . A factor 
1.36 is also included here to account for helium and heavier el¬ 
ements consistent with the SL-sample. For the radial distribution 
of the molecular gas, Eh 2 , Wolfire et al. (2003) apply a Gaussian 
distribution in the galactocentric range between 3 to 6.97 kpc with 
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a peak at 4.85 kpc and a FWHM value of 4.42 kpc. Beyond 6.97 
kpc they use a radial exponential for E h 2 with a scale length of 
2.89 kpc. Same as the SL-sample, our total ISM gas mass column 
density in the Milky Way disk is then E g = Ear + E h 2 ■ 

2.2 Stellar mass column density profiles 

2.2.1 The SL-sample 

Column densities of stellar mass of all galaxies of the SL-sample 
were obtained by surface photometry of infrared images observed 
by the Wide-held Infrared Survey Explorer (WISE; Wright et al. 
2010) in the W1 band at 3.4 pm, the SIRTF Nearby Galaxies 
Survey survey (SINGS; Kennicutt et al. 2003) with Spitzer/IRAC 
at 3.6 pm, and the 2MASS survey (Jarret et al. 2003) in the K s - 
band at 2.2 pm. A procedure similar as described in Leroy et al. 
(2008) is applied. We measure the median intensities in 10 arc- 
sec wide tilted rings at different galactocentric distances. We find a 
tight relation between the WISE/W1 and SPITZER/IRAC intensi¬ 
ties, I 3 .G = I.O 3 / 3 . 4 , and the shape of the radial surface bright¬ 
ness distributions typically agrees very well. The /C s -band im¬ 
ages are not deep enough to reach the outer rings, as opposed to 
the 3.4 and 3.6 pm observations, but for the inner and brighter 
parts of the galaxies where the data overlap, we find a relation 
Ik, = I. 8 / 3 . 4 . We use this relation and the median deprojected 
intensities / 3 . 4 (f) together with a fixed K-band mass-to-light ratio, 
Tf = 0.5 Mq/ L 0 , k (see Leroy et al. 2008), to convert the surface 
brightness profiles to mass column density profiles E*(r). 

In most cases the central parts of the mass column density are 
contaminated by the contribution of galactic bulges, which have a 
different star formation history and chemical evolution than disks. 
Since our goal is to investigate the chemical evolution of the ISM 
and the young stellar population in the disks of these galaxies, we 
subtract the bulge contribution by a applying a bulge-disk decom¬ 
position algorithm. For this purpose, the bulge contribution is rep¬ 
resented by a general Sersic profile 

Etmf ge{r) = Ef, exp[- 6 „((r /r e ) 1/ " - 1 )]. ( 1 ) 

The disk contribution is represented by either a single exponential 
profile or two exponential profiles broken into inner and outer parts 
so that the total mass column density, Et, is given by 

E t (r) = T,buige(r) + E d ,i exp (-r/hi) at r < n, rea k, (2) 

and 

Et(r) = E bui ge (r) + E d , 0 exp (—r/h„) at r > Vbreak■ (3) 

Note that we measure the galactocentric radius r, the effective 
radius r e , the break radius threat, and the the disk scale lengths 
hi, h 0 in units of the iso-photal radius /? 25 - The mass column den¬ 
sities Ej,, E dp, E d,o are given in Mq pc” 2 , and b n is calculated 
as a function of the bulge shape parameter n through the condi¬ 
tion T( 2 n) = 27(2 n, b n ), where T and 7 are the complete and 
incomplete gamma functions, respectively (Graham 2001). All the 
parameters are determined through fitting the observed radial col¬ 
umn densities using the above equations and are summarised in 
Table 1. 

2.2.2 The Milky Way 

The radial stellar mass column density profile of the Milky Way 
disk has been studied recently by Bovy & Rix (2013). From their 
dynamical measurement they derive a surface column density of 


the stellar disk of E*(Ro=8 kpc) = 38 Mq pc” 2 , and show that 
an exponential radial distribution of E* around this galactocentric 
distance with a scale length of 2.1 kpc is a good fit to their ob¬ 
servations. We adopt this fit for the stellar mass in the Milky Way 
disk. 


2.3 Metallicities and metallicity gradients 

2.3.1 The SL-sample 

The standard method to determine metallicities and metallicity gra¬ 
dients of the ISM or the young stellar population in the disks 
of star forming spiral galaxies is the analysis of strong emission 
lines from Hll regions using empirical calibrations which turn 
flux ratios of particular emission lines into oxygen abundances. 
The calibrations of these “strong line methods” are subject to 
large systematic errors and are heavily debated (see, for instance 
Kewley & Ellison 2008 or Bresolin et al. 2009a). An alternative 
method is the quantitative analysis of absorption line spectra of 
individual blue supergiant stars (BSGs) which has been carefully 
tested in the Milky Way and and has now been applied to galax¬ 
ies as distant as 7 Mpc (Kudritzki et al. 2008, 2012, 2013, 2014). 
While this method is much less affected by systematic uncertain¬ 
ties, it has so far been applied to only a handful of galaxies be¬ 
yond the Local Group. Thus, for the investigation of our SL-sample 
we decide to use the homogeneous and comprehensive study by 
Pilyugin et al. (2014), who applied a new carefully tested strong¬ 
line calibration (Pilyugin et al. 2012) on observations of Hll re¬ 
gions of a sample of 130 late type galaxies providing central oxy¬ 
gen abundances and abundance gradients. Since linear radial gra¬ 
dients depend on the distances to the galaxies and in many cases 
the distances used by Schruba et al. (2011) and Leroy et al. (2008), 
which are also what we adopted, are different from the distances 
used by Pilyugin et al. (2014), we have scaled the metallicity gra¬ 
dients by Pilyugin et al. (2014) to the distances by Schruba et al. 
and Leroy et al. 

In order to assess the zero point of the strong line cali¬ 
bration adopted by Pilyugin et al. (2014), we compare the Hll 
metallicities with BSG metallicities for the galaxies NGC 300 
(Kudritzki et al. 2008), M33 (U et al. 2009), M81 (Kudritzki et al. 
2012), and NGC 3621 (Kudritzki et al. 2014). We proceed with 
the comparison in the following way. The BSG metallicities rep¬ 
resent a full range of heavy elements including the iron group 
and oxygen, and are measured relative to the sun on a logarith¬ 
mic scale with a metallicity defined as [Z\ = log (Z/Zq). Since 
the BSG metallicities including the contribution from oxygen are 
measured relative to the sun, we also relate the Hll oxygen abun¬ 
dances, defined as (O/H) = 12 + log[N(0)/N(H)], to the sun 
by introducing the difference A (O/H) = (O/H) — ( 0/H)q, 
where (O/H)® = 8.69 (Asplund et al. 2009). We then calculate 
A [Z] = [Z]bsg — ZA(O///) at three galactocentric distances 
r //? 25 = 0.3, 0.6, 0.9 and obtain A[Z] = 0.13,0.08,0.23,0.15 
dex for NGC 300, M33, M81, NGC 3621, respectively, as the av¬ 
erage over the galactocentric distances in each galaxy. The average 
zero point difference is (A[ZJ) = 0.15 ± 0.06 dex. For NGC 3621 
we have only r/R 25 = 0.65, 0.93 available for the BSGs. For the 
distance to these galaxies and /?25 we use the values given in the 
BSG papers. We apply this zero point shift to the Pilyugin et al. 
central galactic metallicities. We will also investigate the effect of 
the zero point shift on our results later. 

The radial dependence of the logarithmic oxygen abundance 
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Table 1. The SL galaxy sample 


Name 

NGC 

D 

[Mpc] 

R-25 

[kpc] 

(O/H) 0 

[dex] 

d( O/H) 
dr 

[dex kpc -1 ] 

logEf, 
[M 0 pc" 2 ] 

n 

r e 

[R 25 ] 

logSd.i 
[M 0 pc 2 ] 

hi 

[K 25 ] 

log£ d , 0 
[Mo pc" 2 ] 

ho 

[#25 ] 

Inbreak 

[-R 25 ] 

0628 

7.3 

10.40 

8.93±0.01 

0.051 ±0.0024 

2.449 

0.741 

0.044 

2.596 

0.206 

- 

- 

- 

0925 

9.2 

14.20 

8.63±0.02 

0.034±0.0026 

1.577 

0.665 

0.087 

1.796 

0.282 

- 

- 

- 

2403 

3.2 

7.30 

8.63±0.02 

0.051 ±0.0042 

1.770 

0.419 

0.063 

2.383 

0.228 

- 

- 

- 

2903 

8.9 

15.30 

8.97±0.03 

0.032±0.0044 

2.964 

0.817 

0.030 

2.901 

0.154 

- 

- 

- 

3184 

11.1 

11.90 

8.81 ±0.02 

0.015±0.0031 

1.597 

2.114 

0.169 

2.092 

0.645 

3.085 

0.166 

0.477 

3198 

13.8 

13.00 

8.75±0.04 

0.025±0.0044 

2.087 

0.157 

0.045 

2.275 

0.249 

- 

- 

- 

3351 

10.1 

10.60 

8.97 ±0.01 

0.019±0.0028 

3.131 

0.674 

0.052 

2.660 

0.240 

- 

- 

- 

3521 

10.7 

12.90 

8.98±0.04 

0.039±0.0056 

2.522 

2.705 

0.082 

2.921 

0.207 

- 

- 

- 

3938 

12.2 

6.28 

8.94±0.05 

0.089±0.0097 

0.727 

0.272 

0.189 

2.749 

0.295 

- 

- 

- 

4254 

20.0 

14.60 

8.69±0.03 

0.004±0.0026 

2.626 

0.711 

0.076 

3.045 

0.210 

- 

- 

- 

4321 

14.3 

12.52 

8.89±0.03 

0.017±0.0031 

3.079 

0.500 

0.049 

2.661 

0.331 

- 

- 

- 

4559 

11.6 

17.68 

8.68±0.03 

0.021 ±0.0031 

1.361 

1.236 

0.089 

2.313 

0.166 

2.120 

0.196 

0.473 

4625 

9.5 

1.91 

8.73±0.01 

0.041±0.0146 

- 

- 

- 

2.564 

0.347 

- 

- 

- 

4736 

4.7 

5.30 

8.72±0.03 

0.009±0.0222 

3.586 

1.843 

0.083 

2.830 

0.270 

- 

- 

- 

5055 

10.1 

17.40 

9.02±0.02 

0.026±0.0017 

2.648 

1.110 

0.081 

2.814 

0.199 

- 

- 

- 

5194 

8.0 

9.00 

9.03±0.03 

0.025 ±0.0041 

3.132 

0.748 

0.075 

2.811 

0.414 

3.624 

0.186 

0.701 

5457 

7.4 

25.81 

8.86±0.01 

0.025 ±0.0009 

2.402 

1.697 

0.021 

2.487 

0.176 

- 

- 

- 

6946 

5.9 

9.80 

8.87±0.06 

0.034±0.0101 

2.855 

0.584 

0.040 

2.801 

0.275 

- 

- 

- 

7331 

14.1 

19.60 

8.82±0.08 

0.011 ±0.0090 

3.003 

1.460 

0.026 

3.086 

0.136 

2.036 

0.366 

0.490 


is then given by 

(0/H)(r) = 12 + log Hg) = (0/H)o - (4) 

The values of the zero points (0/H)o and linear gradients 
for each galaxy are given in Table 1. 

2.3.2 The Milky Way 

The investigation of the radial metallicity distribution of the ISM 
or the young stellar population in the Milky Way disk requires 
accurate knowledge of the distances of the individual objects in¬ 
vestigated. In this regard, Cepheid stars are by far the best trac¬ 
ers of metallicity as a function of galactocentric distance because 
their distances can be inferred from the application of period lumi¬ 
nosity relationship. Very recently, Genovali et al. (2014) published 
a comprehensive spectroscopic study of 450 Cepheids and deter¬ 
mined the metallicity and metallicity gradient in the same range 
of galactocentric distances for which we have information about 
the gas and stellar mass column density distributions. They ob¬ 
tain [Z] = (0.57 ± 0.02) — (0.06 ± 0.002)r, where r is in kpc. 
We note that Genovali et al. (2014) actually measured the spatial 
distribution of Fe/H but we use this as a proxy for metallicity. 
The metallicity gradient of 0.06±0.002 dex kpc -1 is very well 
constrained by this study, while the zero point could still be sys¬ 
tematically affected because non-LTE effects were neglected. On 
the other hand, Nieva & Przybilla (2012) have provided an ac¬ 
curate determination of cosmic abundance standards by a care¬ 
ful non-LTE investigation of B-Stars in the solar neighbourhood. 
They obtain a total metallicity mass fraction of Z = 0.014 and an 
iron abundunce of 12 + log N(Fe)/N(H) = 7.52, which both are 
very similar to the sun and the protosolar nebula (Asplund et al. 
2009). We use this result to apply a correction to the zero point 
of the Genovali et al. relation by —0.09 dex, resulting in a rela¬ 
tion for the Milky Way disk metallicity of the young stellar pop¬ 
ulation, [Z] = (0.48 ± 0.02) - (0.06 ± 0.002)r, which yields 
[Z] = [Z] B -stars = 0 at the galactocentric distance of r = 8 kpc. 


3 THE CHEMICAL EVOLUTION MODEL 

There is a large variety of model approaches to describe the cir¬ 
culation of matter from the birth of stars in the ISM, through 
their life with stellar winds and mass-loss to their late phases 
as planetary nebulae or supernovae ranging from the simple 
closed-box model (Searle & Sargent 1972) to models with ra¬ 
dial gas flows (Edmunds & Greenhow 1995) or radial stellar mi¬ 
gration (Sellwood & Binney 2002; Schoenrich & Binney 2009a,b; 
Minchev et al. 2013) and models with a chemo-dynamical descrip¬ 
tion (Burkert et al. 1992; see also Matteucci 2012 or Pagel 2009 
for an overview). Since the goal of this study is not to investi¬ 
gate star formation history or the chemical enrichment of the older 
population but solely to focus on the ISM and the young popula¬ 
tion, we will use a simple analytical model. The model includes 
the effects of galactic winds and accretion and allows us to ob¬ 
tain constraints on galactic mass-loss and accretion mass-gain in¬ 
ferred from the ISM and young stars. In the following, we de¬ 
rive the model. We stress that this model is not new and a spe¬ 
cial case of the more general models introduced by Recchi et al. 
(2008). Similar models were also used to study the global metallic- 
ities of galaxies (Spitoni et al. 2010, Dayal et al. 2013, Lilly et al. 
2013, Pipino et al. 2014, Yabe et al. 2015, Zahid et al. 2014). We 
also refer to the pioneering work by Edmunds (1990). 

We consider a volume element in a galactic disk with stellar 
mass M* and gas mass M g and start from the mass balance equa¬ 
tion 

dM g = — dM t — ( Mloss — Maccr)dt. (5) 

The first term on the right-hand side is the gas consumed by star 
formation. We note that M* 1 is the presently observed stellar mass 
and that its increase through the formation of new stars with a star 
formation rate t)i(t) is dM * = (1 — R)ipdt, where R is the fraction 
of stellar mass returned to the ISM through stellar winds and other 

1 Note that the M* here is equivalent to M* 0 in Ho et al. (2015), where 
the subscript “o" denotes the “observed” stellar mass. 
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processes. Note also that we make use of the instantaneous recy¬ 
cling approximation assuming that this return happens on a short 
time scale compared with galactic evolution. The second term de¬ 
scribes the mass-loss and mass-gain through galactic winds and ac¬ 
cretion, respectively; in this formulation both Mi oss and M accr are 
required to be positive or zero. 

Expressing mass-gain and mass-loss in units of the star forma¬ 
tion rate by introducing the mass loading factor 


V = 


Mio 


i/> 


and the mass accretion factor 


A = 


M a . 


and with the definition 


(6) 

(V) 


_ t^tloss A facer t) A 

Q = (1 - R)ip = 1 - R’ 

one can rewrite Equation 5 as 


( 8 ) 


dM g = —(1 + a)dM,. 


(9) 


With the definition of metallicity Z = Mz/M g (Mz is the 
mass of metals in the ISM), one can calculate the derivative of Z 
with respect to M, 


dZ 

dM* 


1 , dM z _ dMg 

Mg 1 dM* dM* 


Mg 


dM z 

dM, 


-f- Z( 1 + cx) 


( 10 ) 


The first term on the right hand side of Equation 10 can be 
expressed as 


dM z _ yz _ „ _ C 
dM* ~ (1 - R) (1 - R) ’ 


where £ is defined as 


* _ m „ dM OSS ry dM a 

ccr 


= Zu 


Mlos 

Ip 


-Z a 


M a 


tp 


— ZlossV Z a ccrA. 


(ID 


( 12 ) 


The assumption Zi OS3 ~ Z has already been discussed by 
Ho et al. (2015). It seems valid if the energy and momentum of the 
sources generating the outflows (for instance, supernovae or stellar 
winds) is transfered to the neighboring ISM so that a fraction of the 
local ISM is lost in galactic winds, a scenario which is supported by 
observations (Rupke et al. 2013), although counter examples like 
the case of the dwarf galaxy NGC 1569 (Martin et al. 2002) have 
also been detected. A more general model, where Z accr and Zi oss 
are arbitrary, could of course also be applied (see Dalcanton 2007) 
but would introduce additional parameters. 

It is worth pointing out that by definition rj and A must be 
positive or zero 


rj > 0 and A > 0, 


(15) 


because of the ways M; oss and Maccr were defined. However, a 
defined in Equation 8 can be both positive (outflow dominates), 
negative (inflow dominates), or zero (balance between inflows and 
outflows). 

To solve the relation between Z, M* and M g , we further as¬ 
sume that R, yz, y, and A (and therefore a) are all constant as a 
function of time and space. While constant R and yz is a standard 
assumption in many chemical evolution models, the assumption of 
constant mass-loading and mass accretion factors requires an ex¬ 
planation. First of all, this assumption allows to obtain an analyti¬ 
cal solution for metallicity as a function of the ratio of stellar mass 
to gas mass. We note that by this assumption we do not restrict the 
time dependence of star formation or mass infall or loss, we only 
assume that the ratios of wind loss or accretion gain to star for¬ 
mation rate are on average constant. Since very likely the energy 
and momentum driving the outflow is related to star formation, r/ = 
const, is a reasonable zero order approximation. At the same time, 
inflow increases the reservoir of gas capable to form stars and, thus, 
a constant A seems also a good starting approximation. The work 
by Recchi et al. (2008) investigating effects of non-constant A con¬ 
firmed this conclusion. Recent work by Bouche et al. (2010) and 
Lilly et al. (2013) also indicates that mass accretion gain is propor¬ 
tional to the star formation rate. 

Given the relation between dM, and dM g in Equation 9, we 

find 


(1 — R) (1 + tx) 
ZA — yz 



(16) 


There are three terms in the right-hand side of Equation 11 , and 
each term describes a different mechanism of increasing or decreas¬ 
ing the metal mass dMz under the formation of some stars dM,. 
The first term represents the metal production by star formation, yz 
is the nucleosynthetic yield defined as the fraction of metal mass 
per stellar mass which is produced by star formation. The second 
term is the fraction of metals locked up in each generation of stars. 
The last terms describes the gain and loss of metals due to inflows 
and outflows. Here, we implicitly assume that the metals are recy¬ 
cled instantaneously and the gas and metals are well mixed. 

If we assume that the accreted gas is pristine ( Z acC r ~ 0), 
and the gas-loss through galactic winds predominately driven by 
supernovae has the same metallicity as the ISM ( Zi oss ~ Z), we 
can express £ in terms of the mass loading factor rj\ 


c 




Mioss 

7p 


Zr], 


which in combination with Equation 10 and 1 1 yields 


dZ 

dM, 


1 

M g {l-R) 


(yz - ZA). 


(13) 


(14) 


We note that this solution is only valid for a different from —1. 
We will discuss the case a = — 1 separately below. Integrating this 
equation and combining the result with the integral of Equation 9 


M s (0) = M g (t) + (1 + a)M,(t), 

we finally obtain the solution for the metallicity 


m = 


yz 

A 


1 - 


where 


1 + (1 + a) 


A 


M,(t) 
Mg (t) 


(17) 


(18) 


(19) 


(1 — i?)(l + a) ’ 
ri > 0, A > 0, and a ^ —1. 

Equation 18 describes the general case with both inflows and 
outflows. However, there is a restriction in this general case for 
models with very strong inflows where a < — 1 or A > 77 + (1 — . 
In such a situation, there are still solutions possible but only for 
small stellar mass to gas mass ratios limited to At,/M g < —1/(1+ 
a). We also note that models, which start with zero gas mass, M g (0) 
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= 0, and develop with strong inflow (a < —1) have the solution 
of constant metallicity Z(t) = yz /A. Models corresponding to this 
case have been discussed by Bouche et al. (2010) and Lilly et al. 
(2013). 

The case a = — 1 which was excluded before can be easily 
addressed. From Equation 17, M g (t) = const, is obtained in this 
case and A becomes 


A = (1 — R) + y. (20) 

With M g (t) = const, we can then directly integrate Equation 14 
and obtain 

/ A M, (t) \ 

Z(t) = 2jLli- e J . (21) 

In the special case whith only outflows and virtually no in¬ 
flows, i.e. y ^ 0 and A = 0, Equation 18 is no longer valid. How¬ 
ever, for this case one can repeat the derivation from Equation 14 
with A = 0 and show that 




when y ^ 0 and A = 0. (22) 

In the alternative special case with only inflows but no out¬ 
flows (A ^ 0 and r\ = 0), Equations 18 remains valid and the 
reader is also directed to the similar models by Koppen & Edmunds 
(1999). 

In the special case where there are no inflows nor outflows 
Equation 22 becomes: 

^) = T^ ln 

This is the classical closed-box solution (Searle & Sargent 1972; 
Pagel & Patchett 1975). We note that in the limit of 1 

all different cases of the models with inflow or outflow (Equa¬ 
tions 18, 21, 22) approach the closed-box case Z(t) = 

Equations 18, 21, 22, 23 describe the evolution of total met- 
alliticity with the change of the ratio of stellar mass to gas mass. 
While in recent years quantitative stellar spectroscopy has ad¬ 
vanced beyond the Local Group and has been applied to measure 
total metallicity including iron group elements from the popula¬ 
tion of young massive stars in galaxies (Kudritzki et al. 2008, 2012, 
2013, 2014; Gazak et al. 2014), the standard method of metallic¬ 
ity determination in the disks of star forming galaxies is still the 
analysis of ISM Hll region emission lines, which yields the num¬ 
ber fraction of oxygen to hydrogen or, more specifically, the oxy¬ 
gen abundance (O/H) = 12 + log[N(0)/N(H)]. If in the equa¬ 
tions above the metallicity yield yz is replaced by the oxygen yield 
yo, then they can be used to calculate the oxygen mass fraction 
Om, which is then converted to the oxygen number fraction by 
N(0)/N(H) = Om/{ 16X) (X is the hydrogen mass fraction). 

The formulae derived, if applied to the presently observed ra¬ 
dial distribution of the ratio of stellar to gas masses in the disks of 
star forming galaxies, predict the metallicity mass fraction of the 
interstellar medium or the young stellar population. The free pa¬ 
rameters of the model are the metallicity or oxygen yield yz or 
yo, the stellar mass return fraction R, and the mass loading and 
accretion factors y and A. We will empirically calibrate yz and yo 
in the next section using Milky Way observations described in the 


1 + 


M,{t) 


Mg(t) 


when y = 0 and A = 0. (23) 


\ M*(f) 

) M g (t) 



Figure 1 . Oxygen abundance as a function of the ratio of stellar to gas mass 
for a characteristic selection of the parameters y and A. Black: closed-box 
model with y = 0, A = 0. Blue: A = 0, y = 0.5 (dashed) and 1.0 (solid). 
Red: A = 0.5, y = 0.0 (dashed), 0.5 (dashed-dotted) and 1.0 (solid). 
Green: A = 1.0, y = 1.0. 


previous sections. The value of R depends crucially on assump¬ 
tions for the initial mass function (IMF) and the strengths of stellar 
winds and mass loss in all phases of stellar evolution and is, thus, a 
rather uncertain number. Leitner & Kravtsov (2011) present calcu¬ 
lations of R as a function of the age of the stellar populations and 
for different IMFs and find a range of approximately 0.15 to 0.5. 
For our study presented here we adopt R = 0.4 throughout. How¬ 
ever, we have also carried out all calculations and fits with a value 
of R = 0.2 in parallel. The resulting changes are small and will be 
discussed in the forthcoming sections. 

The determination of the remaining two parameters y and A 
from the comparison of the observed metallicities and the model 
predictions is the major goal of this paper and will be described 
and discussed in the following sections. At this point, we discuss 
the major properties of the model. Fig. 1 shows how these parame¬ 
ters influence the metallicity as a function of the ratio of stellar to 
gas mass. The closed-box model has the largest metallicity. It does 
not show a saturation in metallicity but continuous with a logarith- 
mical increase as a function of M t /M g . This is a consequence of 
the assumption of a constant yield during the course of chemical 
evolution. Models with A = 0 have a somewhat lower metallicity 
and the metallicity decreases with increasing y. These models do 
not saturate and increase logarithmically. However, models with A 
larger than zero saturate at a metallicity X max — yz I A as long as 
a > — 1. How fast this level of saturation is reached with increasing 
M t /Mg depends on y. Models with very small values of y follow 
the closed-box model closely, but their metallicities suddenly turn 
into the maximum metallicity Z max , whereas models with larger 
y depart early from the closed-box case and then approach Z max 
gradually. 

It is obvious that the systematically different behaviour of the 
curves in Fig. 1 as a function of the parameters y and A opens 
the possibility to determine these parameters for disk galaxies with 
well observed spatially resolved radial profiles of stellar and gas 
mass and metallicity. We demonstrate this in Fig. 2, where we have 
adopted a simple exponential profile for the ratio of stellar to gas 
mass, Mt/Mg = 5 exp (—x/h) with x = r/7?25 and h = 0.5 
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Figure 2. Oxygen abundance as a function of galactocentric radius for a 
galactic disk with an exponential profile of the ratio of stellar to gas mass 
(see text) using the same models as in Fig. 1. 


Figure 3. Oxygen abundance as a function of the ratio of stellar to gas 
mass predicted by the global galaxy model by Ascasibar et al. (2014) and 
described by Equation 24. Blue: T = 4.5, green: T = 8.5, red: T = 11.0. 
The models in grey are the same as in Fig. 1. 


(see Ho et al. 2015 for the choice of this exponential). The closed- 
box model produces the highest metallicity and the steepest gradi¬ 
ent. Models with no inflow are slightly shallower and have lower 
metallicities which depend on the strength of the adopted outflows. 
Models with inflow produce even flatter gradients with increasing 
A. 

In two recent publications, Zahidetal. (2014) and 
Ascasibar et al. (2014) have developed analytical chemical 
evolution models which relate the averaged metallicities of 
galaxies as observed in the large galaxy surveys to the ratio of total 
galactic stellar mass to total gas mass. These models reproduce the 
observed global MZR of galaxies very well. Zahid et al. (2014) 
use the observed scaling relation that the total ISM hydrogen mass 
in the disks of star forming galaxies relates to the total stellar 
disk mass by M g oc with an observed exponent g ~ 0.5. 
However, in the spatially resolved case this relationship does not 
hold and, thus, the approach by Zahidetal. (2014) will require 
some modifications to predict the stratification of spatially resolved 
metallicity in galaxies. 

Ascasibar et al. (2014) apply an analytical time averaging pro¬ 
cedure of the general chemical evolution equations and show that 
their weak wind solution (the case where their ew = 0) 


12 + lo S T77TTT = 9.750 + log( : 


N(H) ' ^Vi + _|_ Ts 

with M g = 1.35 Mh, the ISM hydrogen mass, and 


R M , 
1 — R Mh 


(24) 


(25) 


leads to very good agreement with observed global galaxy oxygen 
abundances as a function of global gas and stellar mass. They ar¬ 
gue that in the way as it was derived this relationship should also 
hold for the spatially resolved case in galaxy disks. Ascasibar et al. 
(2014) suggest a range of their parameter T from 4.5 to 11 and use a 
return fraction R of R = 0.18. Fig. 3 compares the Ascasibar et al. 
(2014) model with three values of T with our models. They all pre¬ 
dict much higher metallicity at low ratios of stellar to gas mass. 
This is a consequence of the higher yield used by Ascasibar et al. 
(2014), which leads to their constant 9.750. Modifying the yield to 


lower values would shift the curves downward. We will test this 
model in the next section. 


4 THE RADIAL METALLICITY DISTRIBUTION OF 
THE MILKY WAY DISK AND CALIBRATION OF 
METALLICITY YIELD 

With the well observed radial distributions of stellar mass, ISM gas 
mass and metallicity in the Milky Way disk as described in Sec¬ 
tion 2, we can now empirically calibrate the metallicity yield yz 
of our chemical evolution model of Section 3. We are aware of 
the many theoretical model calculations of yields (see Zahid et al. 
2012, Yabe et al. 2015, Pagel 2009, or Matteucci 2012 for an 
overview), but with the apparent uncertainties introduced by the 
choices of IMF or assumptions about stellar winds and mass-loss, 
we regard a differential study of galaxies relative to the Milky Way 
with an empirically calibrated yield as more appropriate approach. 

We start with the assumption of the closed-box case described 
by Equation 23. According to Section 2, the stellar mass column 
density of the Milky Way disk is 38 M@ pc -2 and the total gas 
mass of the ISM (Hi and H 2 ) is 8.465 Mq pc~ 2 , so that the ratio of 
the two is 4.49. With a metallicity of the young stellar population 
Zb— stars = 0.014 (see Section 2) and Equation 23, we constrain 
the metallicity yield to yz /{1 — R) = 0.00822. We can use this 
yield and apply Equation 23 on the observed radial profile of the 
ratio of stellar to gas mass (Fig. 4) to predict the metallicity dis¬ 
tribution of the ISM and the young stellar population in the Milky 
Way disk. The result is shown in Fig. 4 and compared to the ob¬ 
served metallicity gradient. The agreement, though not perfect, is 
surprisingly good. 

In the next step, we relax the closed-box assumptions and 
use the full chemical evolution model with galactic winds and ac¬ 
cretion as described in Section 3. We adopt a stellar wind return 
fraction R = 0.4 (see Section 3 for discussion) and a range of 
metallicity yields yz from 0.004 to 0.012. For each yield yz,i 
adopted in this range, we calculate chemical evolution model dis¬ 
tributions [Z] mo d(r, r/j, Afc) for an array of discrete values of r/j 
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MilkyWay MilkyWay 




Figure 4. Left: The radial profile of the observed ratio of stellar to gas mass column density of the Milky Way disk. Right: Observed Milky Way disk metallicity 
[Z] (black) compared to the prediction of the closed-box model (red) and the best fitting chemical evolution model (blue) with r/ = 0.432 and A = 0. 
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Figure 5. isocontours in the fry A)-plane of the observed Milky Way 
disk metallicity distribution [Z] (r) for a yield yz = 0.00671. The Ay 2 
values adopted for the plot are 2.3 (red), 6.17 (blue), 11.8 (black). 

and Afc. We then use the difference between [Z] mo d(r, r/j, Ay) and 
the observed distribution of [Z](r) at hundred radial points dis¬ 
tributed equally between 4.5 to 9.0 kpc to calculate a x 2 -matrix 
X 2 {yz,i, Vj, Ay) and search for the minimum. We find a well de¬ 
termined minimum Xmm at yz = 0.00671, y = 0.432, A = 0. 
The corresponding [Z](r) model distribution is shown in Fig. 4. 
The fit is a clear improvement relative to the closed-box model. 

To assess the accuracy of the yield determination we repeat 
the procedure assuming metallicities at the error margins (see Sec¬ 
tion 2), [Z]+(r) = 0.50 - 0.058r and [Z]-(r) = 0.46 - 0.062r 
(with the uncertainty <r(r) = 0.02 + 0.006r in the \' 2 calcula¬ 
tion). We obtain yz = 0.00776, y = 0.549, A = 0 for [Z] + 
and yz = 0.00572, y = 0.306, A = 0 for [Z)~ and conclude 
that the yield is constrained through our procedure with an accu¬ 
racy of A log yz ~ 0.07 dex. The value of y has the uncertainty of 
Ay « 0.12. While all fits settle at values A = 0, the A\ 2 isocon¬ 
tours at a fixed yield as shown in Fig. 5 indicate an range of AA « 
0.06 to 0.1. 


We also repeat the yield calibration procedure for a stellar 
mass return fraction of R = 0.2 and obtained yz = 0.00895, 
y = 0.576, A = 0 at the best fit solution. The fit is of similar qual¬ 
ity as the one in Fig. 4. We note that with the change from R = 0.4 
to 0.2 the ratios yz/(l — R) and y/(l — R) have remained constant. 
This is a consequence of the fact that the best fits for both R = 0.4 
and 0.2, respectively, have settled on a solution with A = 0 for 
which Equation 22 holds. In the following we will use R = 0.4 
and the correspondingly calibrated yield. 

To apply our chemical evolution model to the observed ISM 
oxygen abundances of the SL-sample of galaxies in the forthcom¬ 
ing section, we need to turn the calibrated metallicity yield yz into 
an oxygen yield yo ■ For this purpose, we use the results of the work 
by Nieva & Przybilla (2012) which we have discussed in Section 2. 
For their sample of B-stars in the solar neighbourhood, they have 
found an oxygen abundance (0/H)B-stars = 8.76. With an aver¬ 
age hydrogen mass fraction X=0.710, this corresponds to an oxy¬ 
gen mass fraction Ob = 0.00654 and an observed ratio of oxygen 
mass fraction to total metallicity mass fraction Ob I Zb = 0.467. 
Thus, in order to obtain with our chemical evolution model the 
same ratio of oxygen mass fraction to total metallicity as observed 
in B-stars in the solar neighbourhood, we need the relation between 
oxygen and metallicity yield as yo = (Ob /Zs)yz = 0A67yz- 
This means that for R = 0.4 we will work with an oxygen yield of 
yo = 0.00313. 

With the good fit obtained in Fig. 4, it is interesting to carry 
out a similar comparison with the chemical evolution model by 
Ascasibar et al. (2014) introduced in Section 3. For this purpose we 
use the relation [Z\ = (O/H) — (0/H)B-stars- The comparison 
shown in Fig. 6 demonstrates that the model which worked well for 
global models of total galactic metallicity using total gas and stellar 
masses does not fit well in the spatially resolved case of the Milky 
Way, because the metallicities saturate too early at still relatively 
low ratios of stellar to gas mass. This is evident from Equations 24 
(see also Fig. 3). The turnover from the linear unsaturated relation 
to maximum metallicity begins at T(R/(1 — R))(M,/M g ) > 1 
or M*/M g > 4.55/T (for R = 0.18, as used by Ascasibar et al. 
2014, larger values of R make things worse). This means that very 
small values of T would be needed to avoid too early saturation. 
Since T in their model is proportional to the ratio of mass-weighted 
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Figure 6. Observed Milky Way disk metallicity [Z] (blue solid) compared 
with the model by Ascasibar et al. (2014) (see Equation 24) for T = 4.5 
(green solid), 8.5 (green dashed) and 11 (green dashed-dotted). 

metallicity of all stars (not only the young population) to the gas- 
phase metallicity, such small values could be obtained by signif¬ 
icant amounts of infall at very early times, which would strongly 
reduce this ratio. However, at the same time, the maximum metal¬ 
licity in their model is given by [Z] ms , x = const. — log T and, 
thus, increases dramatically with decreasing T. To compensate for 
that would require a strong reduction of the yield, as the constant 
in [Z] max changes with A log y z . 


5 ANALYSIS OF THE OBSERVED METALLICITY 
DISTRIBUTION OF THE SL-SAMPLE OF GALAXIES 

With the calibration of the stellar yields obtained in the last sec¬ 
tion, we can now start to analyse the observed radial metallicity 
distribution of our LS-sample of galaxies. We use the observed ra¬ 
tios of stellar to gas mass column densities as a function of galac- 
tocentric radius and calculate model oxygen abundance profiles 
(0/H) mo d(r)(?)i, Aj) with 0 < r]i < 4 and 0 < Ay < 4. For the 
calculation of (0/H) mo d(r), we generally use Equation 18 except 
in those special cases where Equations 21, 22, 23 apply. The dif¬ 
ference of (0/H) mo d(r) to the observed distribution (0/H) o b s (V) 
at discrete radial galactocentric distance points is then used to cal¬ 
culate a x 2 -matrix X 2 (Vi> Ay). The minimum of this matrix then 
defines the best fit of the observed distribution (0/H) o b s (f) and 
provides an estimate of galactic mass-loss y and accretion mass- 
gain A in units of the star formation rate. 

We begin the analysis with NGC2403, which has smooth, 
well observed profiles of Hi, H 2 and stellar mass, and also a good 
coverage of metallicity measurements from Hll region emission 
lines. The result is shown in Fig. 7 and demonstrates that our sim¬ 
ple chemical evolution model can reproduce the observed metallic¬ 
ity distribution surprisingly well. We obtain well constrained val¬ 
ues of galactic wind mass-loss and accretion gain, 77 = 0.81 and 
A = 0.71. 

The galaxies NGC925 and 3198 are similar cases with sig¬ 
nificant amounts of galactic wind mass loss and accretion gain 
(see Table 2). The fit of the observed radial metallicity distribu¬ 
tion in the range where observations exist is reasonable (for all 



eta 


Figure 11. The location of all galaxies discussed in Fig. 7 to 10 in the 
(r/, A)-plane. The green square gives the position of the Milky Way disk as 
discussed in Section 4. 


galaxies we plot the observed metallicities exactly in the range for 
which Pilyugin et al. 2014 have Hll measurements available). The 
increase in the model-predicted metallicity for NGC 925 beyond 
r/i ?25 > 0.9 is caused by a steep decline of Hi gas mass column 
density. 

Contrary to the three galaxies just discussed, we also en¬ 
counter cases with either very small accretion mass-gain or galac¬ 
tic wind mass-loss. For instance, NGC 5457 can be fit very well 
by a model with A = 0, as shown in Fig. 8 . NGC 628 is a case 
with 77 = 0.17. Three more galaxies with very small values of 77 , 
NGC 3184, 3351, and 4559, are displayed in Fig.9. For three of the 
four cases with low galactic wind mass-loss rate, NGC 628, 3184, 
and 4559, the model fails to describe the observed radial metal¬ 
licity distribution for galactocentric radii r/Ris > 0.8, 0.7, and 
0.6, respectively. This is caused by extended and relatively flat dis¬ 
tributions of the ISM Hi gas. The situation is most extreme for 
NGC 4559. Bresolin et al. (2012) and Kudritzki et al. (2014) have 
discussed other even more extreme cases with outer extended disks 
imbedded in large ISM Hi reservoirs as examples where the sim¬ 
ple models used here become insufficient. Our analytical models 
for these three galaxies start to fail when the stellar mass column 
density becomes smaller than the gas column density. 

The galaxies NGC 5055 and 5194 shown in Fig.10 are very 
close to the closed-box case with very small 77 and A values (we 
note that the model distribution for NGC 5194 shows a too strong 
decline at larger radii similar to the cases discussed before). 

We summarise the results for all galaxies displayed in Fig. 7 
to 10 in Fig. 11. The intriguing result of this figure is the indication 
of the existence of three different groups in our sample: 1 ) galaxies 
with either mostly winds and almost no accretion, 2 ) mostly accre¬ 
tion and weak winds, and 3) galaxies with significant amounts of 
both winds and accretion, where the rates of mass-loss and mass- 
gain are of the same order. We also note that according to our fit the 
Milky Way belongs to the subgroup with a significant amount of 
galactic wind but with only a small accretion rate (see also Fig. 5). 

The fits of the remaining galaxies are shown in Fig. 12 to 14. 
For all these galaxies, except NGC 2903, the fits of the observed 
metallicity distributions are reasonable within the uncertainties of 
the observations and the range of galactocentric distances where 
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Figure 7. Chemical evolution fits of the galaxies NGC2403 (top), 925 (middle), 3198 (bottom). Left: comparison of the observed radial oxygen abundance 
distribution (black with the lcr uncertainty indicated by the dashed lines) with the final model (blue) and the closed-box model (red). Right: A\" isocontours 
in the (rj. A)-plane corresponding to Ax 2 = 2.3 (red), 6.17 (blue) and 11.8 (black). 
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Figure 8 . Chemical evolution fits of the galaxies NGC 5457 (top) and 628 (bottom). The arrangement of the figures is similar to Fig. 7. The observations for 
galactocentric distances shorter than indicated by the dashed vertical line are not included in the x 2 -fit. 


observations are available. But for most of these objects the A^ 2 
isocontour diagrams indicate that the values of rj and A are not very 
well constrained by the fit. We add these objects to the plot of the 
location of galaxies in the (r/, A)-plane in Fig. 15 and find that we 
still can identify three well distinguished groups in our sample, in 
the same way as already discussed above. The galaxy NGC 2903 
has an unusual Hi profile with a steep decline of gas mass from 
the centre to 0.5r/J?25 and then an extended disk with constant 
gas mass. Our models are unable to reproduce chemical profiles for 
such gas distribution properly. 

The observed ISM Hll region metallicities taken from the 
compilation by Pilyugin et al. (2014) have been re-calibrated in 
Section 2 using metallicity determinations obtained by stellar spec¬ 
troscopy. These metallicties are given in Table 1 and have been used 
to obtain the results described so far. To investigate the system¬ 
atic effects of metallicity measurements on our method, we repeat 
the analysis described before with logarithmic oxygen abundances 
shifted by A(0/H)o = -0.05,-0.10, and —0.15, respectively, 
going back in small steps to the original zero point calibration by 
Pilyugin et al. (2014). We note that changing the zero points of the 
observed metallicities is equivalent to changing the yield in the 
chemical evolution model (see Equations 18, 21, 22, 23). Thus, 


the results discussed below can also be seen as a test of the influ¬ 
ence of the uncertainty of the yield. 

The results for the different metallicity zero points are given 
in Table 2. Fig. 16 displays the shifts in the (r),A)-plane caused 
by these systematic changes of the metallicity zero points. The 
changes obtained in this way are relatively large in the sense that 
systematically decreasing the metallicty leads to larger values of rj 
and A, which is easily understood with the help of Fig. 1 and 2. This 
changes the conclusions somewhat. While galaxies with very small 
values of A at A(0/H)o = 0.0 remain low accretion cases, most 
of the weak galactic wind objects are now fit with higher mass-loss 
rates. This increases the number of objects with both wind and ac¬ 
cretion, at least for the largest zero point shift A(0/H)o = —0.15. 

Finally, we test the influence of the stellar mass return factor. 
While we argued in Section 3 that R = 0.4 is the preferable choice, 
R = 0.2 cannot be completely excluded. In consequence, we re¬ 
peat our analysis for R = 0.2 and the corresponding yield (see 
Section 4). The results are shown in Fig. 17. While for R = 0.2 
the ?) and A values for the fits obtained are somewhat larger, the 
figure clearly indicates that none of the conclusions made before 
are affected. 
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Figure 9. Chemical evolution fits of the galaxies NGC3184 (top), 3351 (middle), 4559 (bottom). 


6 DISCUSSION AND SUMMARY 

In the previous sections, we have developed and tested a new 
method to use the radial metallicity distributions of the ISM and 
the young stellar population in connection with the radial profiles 


of stellar mass and ISM gas mass column densities to constrain the 
strengths of galactic wind and accretion processes. The method is 
based on the application of an analytical chemical evolution model, 
which assumes constant ratios of galactic mass-loss and mass-gain 
to the star formation rate and which then predicts radial metallicity 
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Figure 10. Chemical evolution fits of the galaxies NGC 5055 (top), and 5194 (bottom). 
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profiles as a function of the ratio of stellar mass to ISM gas mass. 
Generally, this relatively simple model reproduces the observed 
metallicies and their radial gradients very well. In some cases, how¬ 
ever, the model fails at outer galactocentric distances mostly for 
galaxies with extended disks of high Hi gas column densities, but 
it still remains useful in the inner parts. 

We have applied the model to the observations of the Milky 
Way disk and empirically determined the metallicity yield yz . We 
obtained yz/{ 1 — R) = 0.0112 ± 17% = 0.8 %B-stars, where 
%b— stars = 0.014 is the metallicity mass fraction of young B-stars 
in the solar neighbourhood (Nieva & Przybilla 2012). R is the frac¬ 
tion of stellar mass returned to the ISM and has been adopted as R = 
0.4 in our model calculation (other values of R were also tested but 
the influence of varying R was found to be small). Using a B-star 
mass fraction of oxygen to total metallicy O m /^B-stars = 0.467, 
we determined an oxygen yield of yo = 0.00313. With a galactic 
wind mass-loading factor rj = Mi oss /ip = 0.432 ± 0.12 and a 
very low accretion gain A = M acc7 ./t/> = 0.0±o;oo, the model 
reproduces the observed metallicity gradient very accurately. We 
conclude that based on our model fit the Milky Way disk is charac¬ 
terised by a low accretion rate and a moderate galactic wind mass 
loading factor. We note, however, that this result is in disagreement 
with the estimate of the Milky Way infall rate from ultraviolet ab¬ 


sorption line studies of halo high velocity clouds in the sightline 
of halo stars from which Lehner & Howk (2011) estimate an in¬ 
fall rate of 0.8 to 1.4 Mq yr -1 . With a Milky way star formation 
rate of 1.9 ± 0.4 Mg yr -1 (Chomiuk & Povich 2011; Davies et al. 
2011), this would lead to a value of A between 0.35 to 0.93 much 
larger than our very good fit indicates. However, we also note that 
Chomiuk & Povich (2011) do not rule out star formation rates a 
factor of two higher for alternative choices of the initial mass func¬ 
tion and that the determination of the infall rate depends on a va¬ 
riety of parameters which are not well constrained as discussed by 
Lehner & Howk (2011). 

We also tested the alternative chemical evolution models de¬ 
veloped recently by Ascasibar et al. (2014) to reproduce the ob¬ 
served global metallicities of galaxies by the use of the global 
galaxy gas and stellar mass. While this model has been successfully 
applied to explain the observed global mass-metallicity relationship 
of galaxies, it fails in the spatially resolved case of the Milky Way 
to reproduce the observed metallicity gradient. 

About « 30% of our sample of galaxies studied showed simi¬ 
larly low values of accretion rates A as the Milky Way disk with si¬ 
multaneously a range of mass loading factors between 0 < y < 1 . 
However, we also find a group with very low mass loading factor 77 
and 0.3 < A < 0.7. This group comprises ss 40% of our sample. 
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Table 2. 7 / and A determination for different metallicity zero point shifts A(0/H)o 


Name 

A(O/H ) 0 

A(O/H ) 0 

A(O/H ) 0 

A(O/H ) 0 

A(O/H ) 0 

A(O/H ) 0 

A(O/H ) 0 

A(O/H ) 0 


0.0 

0.0 

-0.05 

-0.05 

- 0.10 

- 0.10 

-0.15 

-0.15 

NGC 

V 

A 

V 

A 

V 

A 

V 

A 

0628 

0.17 

0.26 

0.39 

0.25 

0.70 

0.20 

1.07 

0.13 

0925 

0.67 

0.54 

1.10 

0.49 

1.65 

0.39 

2.38 

0.21 

2403 

0.81 

0.71 

1.20 

0.75 

1.64 

0.80 

2.19 

0.83 

2903 

0.00 

0.29 

0.00 

0.38 

0.00 

0.47 

0.00 

0.55 

3184 

0.00 

0.41 

0.07 

0.58 

0.25 

0.52 

0.46 

0.56 

3198 

0.27 

0.52 

0.50 

0.56 

0.76 

0.61 

1.10 

0.65 

3351 

0.00 

0.29 

0.08 

0.33 

0.23 

0.36 

0.42 

0.39 

3521 

0.15 

0.30 

0.40 

0.26 

0.74 

0.18 

1.18 

0.05 

3938 

1.04 

0.00 

1.36 

0.00 

1.74 

0.00 

2.18 

0.00 

4254 

0.00 

0.48 

0.00 

0.64 

0.00 

0.75 

0.06 

0.86 

4321 

0.34 

0.00 

0.55 

0.00 

0.78 

0.00 

1.06 

0.00 

4559 

0.01 

0.66 

0.17 

0.74 

0.41 

0.82 

0.67 

0.91 

4625 

0.40 

0.50 

0.57 

0.56 

0.79 

0.62 

1.04 

0.69 

4736 

0.00 

0.53 

0.00 

0.60 

0.04 

0.68 

0.12 

0.76 

5055 

0.09 

0.05 

0.30 

0.00 

0.51 

0.00 

0.76 

0.00 

5194 

0.00 

0.00 

0.00 

0.04 

0.00 

0.23 

0.00 

0.38 

5457 

0.74 

0.00 

1.00 

0.00 

1.31 

0.00 

1.67 

0.00 

6946 

0.57 

0.00 

0.82 

0.00 

1.10 

0.00 

1.44 

0.00 

7331 

0.00 

0.44 

0.00 

0.53 

0.11 

0.60 

0.27 

0.67 



Figure 15. The location of all galaxies investigated in Fig. 7 to 10 in the 
(r/, A)-plane. The galaxies of Fig. 7 to 10 which have more accurate values 
of 77 and A are plotted as red filled squares and the remaining galaxies as 
blue open squares. The filled green circle gives the position of the Milky 
Way disk as discussed in Section 4. 


The remaining galaxies have roughly equal rates of accretion and 
mass-loss in a range 0.2 < A « 77 < 0.8. 

The range of galactic mass-loss found in our study is in agree¬ 
ment with spectroscopic studies searching for outflows in the local 
universe (e.g., Rupke et al. 2005 ; Veilleux et al. 2005; Bouche et al. 
2012; Zahid et al. 2012; Yabe et al. 2015). Bigiel et al. (2014) iden¬ 
tify clear inflow signatures form their study of Hi velocity fields for 
NGC 2403, 3198, 2903, 7331, for which we also determine signifi¬ 
cant rates of mass infall (see also de Blok et al. 2014). Bigiel et al. 
(2014) also detect infall for NGC 6946. We note that while our fit 
for NGC 6946 finds a solution with A = 0 as the best value, the 
shape of the isocontours in Fig. 14 does not rule out a value of 


Figure 16. The shift of the location of the galaxies of Fig. 7 to 10 in the 
( 77 , A)-plane with different assumptions for the zero point of the metallicity 
determination from Hll regions. Red filled circles refer to the metallicities 
given in Table 1, which were also used for the previous plots. Blue open 
circles are the results when a zero point shift of A(0/H) = —0.15 is 
applied to the logarithmic oxygen abundances. Blue filled squares represent 
shifts of A(0/H) = —0.05 and —0.10, respectively. 


A as high as 0.4. For NGC 925, Bigiel et al. (2014) find no signs 
of inflows but instead significant outflows. With 77 = 0.67 and 
A = 0.54, we detect in NGC 925 both outflows and inflows, but 
the inflow rate within 3 a can still be consistent with zero (Fig. 7). 

The division into three groups with either mostly winds and 
weak accretion or mostly accretion and weak winds or winds and 
accretion roughly equal is an intriguing result. In particular, the 
existence of a significant group of objects disconnected from the 
cosmic web with A « 0 is somewhat surprising, although cos¬ 
mological simulations do not rule out such cases in the local uni- 
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Figure 12. Chemical evolution fits of the galaxies NGC3521 (top), 3938 (middle), 4254 (bottom). 


verse (Keres et al. 2005; Bouche et al. 2010; see Yabe et al. 2015). 
Our results do not confirm the observational finding by Yabe et al. 
(2015) based on an investigation of global metallicities, masses and 
star formation rates that galaxies evolve along the “equilibrium re¬ 
lationship” A = rj + (1-R) by Dave et al. (2011a). It is interesting 


that Yabe et al. (2015) applied the same chemical evolution model 
as this study; however, an important difference is that in our investi¬ 
gation the spatially resolved information about metallicity, gas and 
stellar masses is used and the gas mass profiles are obtained directly 
from radio observations, whereas in Yabe et al. (2015) the infor- 
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Figure 13. Chemical evolution fits of the galaxies NGC4321 (top), 4625 (middle), 4736 (bottom). 


mation about gas masses is obtained indirectly from the Kennicutt- 
Schmidt law and Ha observations. In addition, Yabe et al. (2015) 
do not disentangle the effects of galactic bulges by a bulge decom¬ 
position of the stellar mass profiles. 

We find that the members of the three galaxy groups have 


distinct properties of their gas disks. NGC2403, 925, 3198, 4625 
which comprise the group with relatively large and roughly equal 
r/ and A values in the diagonal of Fig. 15 are all characterised by 
relatively metal-poor Hi dominated gas disks with much lower H 2 
column densities and an almost flat radial profile of the total gas 
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Figure 14. Chemical evolution fits of the galaxies NGC6946 (top), 7331 (middle), 2903 (bottom). 


mass column densities (AlogE g < 0.3 dex R^jj 1 ; see figures in 
Schruba et al. 2011). Of the remaining galaxies only NGC4559 
has similar gas disk properties. We note that NGC4559 has also 
a high A value (0.66) but a very small value of 77 . The gas disks 
of all other galaxies are dominated by molecular hydrogen with 


column densities larger than the one of Hi in the inner disk re¬ 
gion (r < 0.5 to 0.7 R 25 ) and show a clear exponential decline 
(A log Eg > 0.6 dex of their total gas mass column den¬ 

sities over the galactocentric radius range used for our metallicity 
fit. It seems that among these galaxies the ones with low accretion 
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eto 


Figure 17. Same as Fig. 16 but for for the stellar mass return fraction R = 
0.2 (see text). 


rates have on average a steeper exponential decline. For the low 
A cases NGC5194, 5055, 4321, 6946, 5457 we find AlogEg = 
1.4,1.21.3,1.2,1.1 dex R^, respectively (NGC3938, another 
low accretion case but with a very high tj = 1.04, is an excep¬ 
tion with A log E g = 0.7 dex R^s)- On the other hand, most 
of the galaxies for which we determine low outflow rates, have 
mostly flatter profiles. NGC628, 3184, 3351, 3521, 7331 show 
A log Eg = 0.6 to 0.7 dex R^j. The exceptions are NGC4736 
and 2903 with a steeper decline. NGC 4254 also shows a high value 
of A log Eg = 1.3, but the isocontour in Fig. 12 does not rule out 
this galaxy as a potential low accretion case. 

For the Milky Way, we have also determined a relatively low 
accretion rate (Figs. 5 and 15), however, its ISM gas properties are 
different from the rest of the low accretion group. In the range of 
galactocentric distance from 4 to 9 kpc, which we used for the fit 
in Fig. 4, the Milky Way gas disk is dominated by atomic hydrogen 
and has a relatively flat radial profile. The stellar mass column den¬ 
sity is much higher than the gas mass column density. Comparing 
the Milky Way to the group in the upper diagonal of Fig. 15, for 
which the gas disks are also flat and dominated by Hi, the metallic- 
ity of the Milky Way is significantly higher. 

As any other study of this kind using metallicities obtained 
from Hll region emission lines, our results depend on the cal¬ 
ibration of the strong-line diagnostics adopted. In our work, 
we have re-calibrated the Pilyugin et al. (2014) oxygen abun¬ 
dances by applying a shift A(0/H)o = 0.15 dex based on the 
quantitative spectroscopy of blue supergiants in spiral galaxies 
(Kudritzki et al. 2008, 2012, 2013; U et al. 2009). If we use the 
original Pilyugin et al. (2014) calibration, we obtain larger values 
for t] but the qualitative conclusions remain the same as before. In 
particular, the fraction of galaxies with very small accretion rates 
A remains unchanged. For future work, we regard it as crucial to 
improve the calibration of the Hll region strong line methods based 
on a larger sample of stellar spectroscopy of blue and also red su¬ 
pergiants (Gazak et al. 2014). 

In summary, we conclude that the use of spatially resolved in¬ 
formation of observed metallicity, stellar mass and gas mass is a 
powerful tool to constrain the effects of galactic winds and accre¬ 
tion on the chemical evolution of star forming galaxies. 
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